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1 Introduction 

Semiclassical quantization is very useful to under- 
stand the global behaviour of eigenfunctions and 
energy spectra of quantum systems, since it al- 
lows us to obtain analytic expressions. The lead- 
ing semiclassical approximation (torus quantiza- 
tion) is just the first term of a certain ^-expansion, 
which is called WKB (Maslov and Fedoriuk 1981). 
A systematic study of the accuracy of semiclassical 
approximation is very important, especially in the 
context of quantum chaos (Casati and Chirikov 
1995, Gutzwiller 1990). Since this is a difficult 
task, it has been attempted for simple systems, 
where in a few cases even exact solutions may be 
worked out (Bender, Olaussen and Wang 1977, 
Voros 1993, Robnik and Salasnich 1997a,b, Salas- 
nich and Sattin 1997). 

In this paper we analyze the energy splitting of 
doublets in a generic one-dimensional double-well 
potential. By using the WKB quantization we de- 
duce an analytical formula for the energy splitting 
which is the usual Landau (1997) formula with 
additional quantum corrections. The splitting for- 



mula can be formally written as 
AE = A exp 



(1) 



where S is the usual classical action inside the clas- 
sically forbidden region (between the two turning 
points) and A is called the tunneling amplitude, 
which can be written as a polynomial expression, 
expanded in powers of Ti. This formula is based 
on a linear approximation of the potential near 
the turning points. First we introduce the basic 
definitions, then we derive in detail the splitting 
formula by using the semiclassical WKB expan- 
sion and finally we study its validity for the dou- 
ble square well potential, the inverted harmonic 
oscillator and the quartic potential. Our present 
work is essentially a semiclassical expansion the- 
ory of the tunneling amplitude A of equation (1). 
We shall demonstrate that it is indeed a signif- 
icant improvement of the Landau approximation 
(1997). Also, we shall mention some potential ap- 
plications. 
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2 Basic formalism 

Let us consider a one-dimensional system with 
Hamiltonian 



(2) 



where V(— x) = V(x) is a symmetric double-well 
potential. The stationary Schrodinger equation of 
the system reads 

t^2 j2 

H^(x) = {- — ^ + V(x))ij(x)=E^{x). (3) 

The Sturm-Liouville theorem (see, for example, 
Courant and Hilbert 1968) ensures that for one- 
dimensional systems there are no degeneracies in 
the spectrum. Let ipi and -02 be two exact eigen- 
functions of the Schrodinger equation 



Hipi = Eiipi and Hip2 = ^2^2 



(4) 



such that ipi(-x) = ipi(x) and ^(—x) = —ip2(x) 
and Ei ~ E 2 . To calculate the splitting AE = 
E2 — Ei , we multiply the first equation by ip2 and 
the second by ipi and then we subtract the two 
resulting equations. By integrating from to 00 
we find 



AE 



h 2 Vi (0)^2(0)- ^1 (0)^(0) 



2m 



Io° ipi(x)ip 2 (x) dx 



(5) 



We write the eigenfunctions ipi and ip2 in terms of 
the right-localized function 

Mx) = ^(Mx) + M^))- (6) 

It is easy to show that Eq =< ipo\H\ipo >= ^{E\ + 
E2). Then, with the approximation J* °° ip^dx ~ 1, 
committing an exponentially small error, namely 

roo poo I Y 

if) 1 (x)ip2(x)dx = J ij%(x)dx- - » - , (7) 



we get 



2ti 



AE = VoWo(O) , 



m 



(8) 



which is an almost exact starting formula to calcu- 
late the energy splitting, since the error committed 



in approximation (7) is exponentially small. One 
should observe that this quantity is always posi- 
tive, because the tail of the right localized eigen- 
function ipo(x) at x = has the same sign for V'o(O) 
and its derivative ip'o(0). Another way to see this is 
to realize that due to the Sturm-Liouville theorem 
there are no degeneracies in one-dimensional sys- 
tems, implying that all pairs of almost degenerate 
states, from the ground state up, are grouped by 
odd state above the even state. 

3 Semiclassical method 

To determine the function ipo we perform a 
WKB expansion of the Schrodinger equation. We 
observe that a generic eigenfunction ip of the 
Schrodinger equation can always be written as 



V>(x) = exp(-a(x)) , 



(9) 



where the phase a(x) is a complex function that 
satisfies the Riccati differential equation 



a' 2 {x) + (-)a"{x) = 2m(E - V(x)) . (10) 



The WKB expansion for the phase is given by 



00 ^ 



(11) 



fc=0 



Substituting (11) into (10) and comparing like 
powers of Ti gives the recursion relation (n > 0) 
(see Bender, Olaussen and Wang 1977) 

n 

a'l = 2m(E - V{x)) , ]T a' k a' n _ k + a 1 ' ' n _i = . 

k=0 

(12) 

With the momentum p = ^j2m{E — V(x)) the 
first five orders in the WKB expansion are given 
by 



^0 = P: 



(Jo = 



2p' 
Ap 2 8 p 3 
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4 



0-4 



p'" 3 p'p" 3 p' 3 

J_ /p™ _ gV 13 p" 2 ^ 



16 I p 4 



2 p 5 



+ 



1 /99 p"p 1 297 p 
16 IT p 6 ~ 



'4 N 



(13) 



4 



1 / p w " ir p"V 



32 V P 5 



+ 15 x - 



+ 24 



32 



/// '2 "2 / v 

-111**- - 144^- 



p> 



p' 



+L ( 510^ - 306^ 



32 



p' 



In particular, if we call a and b the two turning 
points corresponding to the energy E, the right 
localized wavefunction tpQ is given by 

Ci 



C 2 

H 7F=F ex P 



cxp 



\p\ dx + (7, 



1 

' i 11 I T~ ~~j~^a,l,l 



\p\ dx + 



even 



1 

r 



for a < x < b (allowed region), where 

00 

CTeven = ^{-l) k H 2k (T 2k (\p{x)\) 
k=l 



(14) 



with a 2k (-\p\) = -&2k(\p\), 



(15) 



and 



a odd = ^(-l) fe n 2fc+ V 2jfc+1 (b(x)|) 
fc=i 

with a 2k+1 (-\p\) = <T2k+i(\p\) ■ 
Instead we get 

C a 

- ran 

h 



(16) 



V>o(x) 



: CXp 



- (^J \p\dx + a even + a odd 



(17) 



for x < a (forbidden region), and also 



ip (x) 



C b 



= cxp 



^ (~ lb dx ~ ^ et,en + a ~ odd ) 

(18) 



for 6 < x (forbidden region), where 



E^2 fe (b»|) 

fc=l 



and 



with a 2k {-\p\) = -°2k 



Oodd = Y. %2k+l ^k+i{\p{x)\) 

k=l 



(19) 



with a 2k+1 ( 



a 2k+1 (\p\), (20) 



where p = y/2m(V(x) — E) and (?' k (x) = a' k (p(x)). 
In evaluating the integrals o k using (13) we cannot 
integrate naively on the real axis, because such 
integrals are divergent, but must take a partial 
derivative w.r.t. the energy of certain fundamental 
complex contour integral. See (Bender el al 1977) 
and section 6. 

We observe that it is possible to write C a , Cb, 
C\ and C 2 in terms of a unique parameter C by 
imposing the uniqueness of the wavefunction ipo 
at the turning points. Following Landau (1997) 
and Merzbacher (1970) we suppose that near the 
turning point x = a it is possible to approximate 
the potential locally linearly by writing 



E - V(x) = F (x - a) , 



(21) 



with Fo > 0. In this case the connections at the 
turning point imply that 



C a = C b = C , d = Ce*4 , C 2 = Ce 



(22) 



and the right localized function i]}q(x) can be writ- 
ten for \x\ < a as 



tpo(x) 



C 



= exp 



w\ 



dx + a even + a dd 



(23) 
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In this way the splitting formula, up to the 5th 
order, after taking into account some straightfor- 
ward relations for and its derivatives (see the 
Appendix), becomes 



AE 



2hC 2 



m 



1 + fi,2<=r 2 (0) + h 2 2dj(0) 



+h 3 (W(0) + ^ 2 3 (0)) 



+ft 4 (-^(0)+4<? 2 (0)a 4 (0) 



(24) 



x exp 



2 

h Jo 



\p\ dx 



The equation (24) can be written in a compact 
form, certainly up to the 5th order, and probably 
also generally to all orders, as follows 



. „ 2hC 2 
AE = exp 



m 



dx + a even (0) 



M= f a' even (\p{x)\)dx. 

J a 



(25) 



To determine C we impose the normalization con- 
dition 

/•oo 

1= / \Mx)\ 2 dx, (26) 
Jo 

from which we get 

2C 2 [ b -exp(2^)dx = l, (27) 



and an expression for C 2 



C 2 



+2h A 







( - dx 




Ja P 





i + 2h 



f b dx 
2 -la P X 



fa f dx 



Jaldx 



rb "j; 
Ja 

Jaidx 



2 -b o-2+(T5 



dx 



(28) 



which is valid up to the 4th order in the T% power 
series of C 2 . The final formula is given by 



AE = — 

m 



' ,6 i 

J a P 



-i -1 



dx 



f b £a dx ' 

x {1 + 2ha 2 (0) + 2h 2 ( a 2 (0) + J -^j— 

la p dx 

f b ^ dx' 

+2h* \ cr 4 (0) + ^(0) + 2^ 2 (0)^f— 

Ja p d x 



+2h A f-a 4 (0) + 2a 2 (0)a 4 (0) 



+2h 4 2or|(0) 



fn — dx 

Ja v 



Hldx 



(29) 



+2h A 



x exp 



fa f dx 



fa'pdx 



J a 



0-3+0-5 



dx 



\p\ dx 



Z\dx 



This formula is the usual Landau (1997) formula 
for the energy splitting (1st order in fi for the tun- 
neling amplitude) with additional quantum cor- 
rections (up to the 5th order in h for the tunnel- 
ing amplitude). We note that higher-order WKB 
corrections quickly increase in complexity (Rob- 
nik and Salasnich 1997a,b) but, in principle, they 
can be calculated from the equation (12). It is 
important to stress that our splitting formula is 
good if the potential is sufficiently smooth so that 
the linear approximation is valid near the turning 
points. 

The same splitting formula (25) and (29) can be 
derived using the semiclassical scattering formal- 
ism, for example as expounded by Iyer and Will 
(1987) and by Will and Guinn (1988). However, 
the main result (25) and (29) of this paper cannot 
be obtained by simple substitution or reinterpre- 
tation of their results. They use the scattering ap- 
proach and treat the scattering problem (asymp- 
totical free motion), calculating the transmission 
coefficients for the tunneling penetration through 
a potential barrier, and specifically they treat the 
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behaviour near to the top of the potential bar- 
rier. Also, their approach is doubly perturbative, 
namely they make the power expansion of the po- 
tential around the top of the barrier whose leading 
term is of course the inverted harmonic oscillator 
plus power terms in the series expansion, and they 
do at the same time the semiclassical expansion 
in terms of the powers of the Planck constant fi. 
They offer a formalism (an algorithm) how to cal- 
culate the requested quantities (transmission coef- 
ficients) to all orders, however they solve the rel- 
evant equation only up to the fourth order. So, 
their result cannot be easily mapped (by substitu- 
tions and other simple operations) onto our prob- 
lem and our solution. Indeed, to get the result 
(25) and (29) using the scattering approach it is 
necessary to go back to the very first step in their 
formalism. 

We have done this and confirmed, as mentioned 
above, that the result is the same. To this end the 
solution is written down in the form of the semi- 
classical ansatz in each of five regions separated 
by the turning points. At the turning points we 
do not request the condition of the continuity of 
the wavefunction and its derivative, but use the so- 
called Kramers correspondence rules instead (they 
determine the coefficients of the ansatz in such 
a way, that exponentially increasing solution in 
the classically forbidden regions does not occur). 
Thus, the asymptotic boundary conditions for the 
scattering problem are automatically taken into 
account: no propagation to the left or to the right 
of the classically allowed region. 

Using this ansatz we have six unknown coeffi- 
cients plus the eigenenergy that we seek, and we 
have six linear equations plus the normalization 
condition, so the problem is well defined. This sys- 
tem of equations can be reduced by simple elimi- 
nation of some of the coefficients to a set of two ho- 
mogeneous linear equations for which the solvabil- 
ity condition is now vanishing of its determinant, 
which must be satisfied precisely at the eigenvalues 
of the energy. Assuming that the pair of almost 
degenerate levels is separated by a small amount 
AE, we can do the Taylor expansion up to the 
first order in AE, neglect the quadratic and higher 



terms, and obtain exactly the equation (25) and 
(29). 

4 Double square well potential 

As the first example, we consider the double 
square well potential. In this case the linear 
approximation of the potential near the turning 
point is not valid. The potential is given by 



V(x) 



for \x\ < a 
for a < \x\ < b 
for \x\ > b 



(30) 



For this potential we have p'(x) = p"{x) = for 
—a < x < a and the corrections to the Landau 
(1997) formula are zero. A naive application of 
the splitting formula gives 



AE 



2hy/E 



'2m(b — a) 



exp 



^^2m(V -E) 



(31) 

This formula is not correct. In fact, by using the 
exactQ wavefunction 



V>o (x) = D exp 



Ti 



2m(V -E)x) (32) 



for < x < a (forbidden region), and 

ipo(x) = A exp ^-V2m£ij +B exp ^ — —y/2mEx^ 

(33) 

for a < x < b (allowed region), and by imposing 
the exact matching and normalization conditions 
(Fliigge 1971) we find 



A 



D 




xexp (^y / 2m(V c 



o 



E) 



(34) 



1 Actually, strictly speaking, this is not exact but nev- 
ertheless the same expression that we get by evaluating 
the leading term for AE by using the implicit trigonomet- 
ric eigenvalue equation (Fliigge 1971, Robnik and Salasnich 
1997, unpublished). 
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e X p(y2m(V -E) + ^VE), 



(35) 



and 



D 2 



2E 



— exp ( - ^^2m(V - E)j . (36) where 



V (b-a) 
Then we obtain: 



AE 



4hE^/2m{V - E) 
mVo(b — a) 



exp[-^y/2m(V -E)). (37) 



This is the exact energy splitting for the dou- 
ble square well potential. It differs by a factor 
4:y/E(V - E)/V from the WKB result based on 
the connection formulae (21-22) which are not jus- 
tified in the present case. 

5 Inverted harmonic oscillator 

In this section we compare the Landau formula 
of the energy splitting with the exact one. We 
consider the inverted harmonic oscillator given by 

Vb(l - for \x\ < a 

V{x) = { for a < \x\ < b (38) 

oo for |x| > b 

We can introduce the following reduced variables 

h 



x T b 
X = — , b = — 
a a 



E 



T7 > n eff 



(39) 

Then it is not difficult to show that the Landau 
formula reads 



AE = 2h eff D 2 exp ( 
where 



7T 



h eff V2 



(l-E 



D 



-1 + VE\ (6 — 1)-| 



1 fl + VE\ 



Ve 



(40) 



(41) 



This formula can be compared with the exact en- 
ergy splitting. Let ip(x) be a quantum state of 
the inverted harmonic oscillator. It can be writ- 
ten in terms of the eigenstates and ifj^ of 
the square well potential 



^(x) = <T(a n ^+b n ^),, (42) 



71=1 



^4+) = -j= cos 



2n - 1 

7T X 



^ ) = Vb 



2b 

n 

TT — X 

b 



The matrix elements of the quantum Hamiltonian 
of the inverted harmonic oscillator read 



J —b 

f m = n: % (^f + & 
4 



= < 



" ) r»(2n-l)» &m 
" (2n-l) 2 7r 2 COS 



2n-l 



b 

2n-l 



m j£ n : 



7T 
7T 

7r 3 I (m—n) 3 (ra+n— l) 3 



26 / cos 



m — n n 
b 



(m—n) 2 



+ 



(m+n-1) 2 



(43) 



and 



*4~2 = j\k ] H^dx 



m = n 



Si (H^) 2 + -2. 
2 U 71 / ^36 



26 2 



sm 



7r 3 (2rt) 3 

2fe „„„ 
"(2n) 2 7r 2 COS 



2n 



7T 



m ^ n 



2b 2 



2b 



(m—n)' 3 



(m—n) 2 



(m+n) 3 



„ r m + n 1 

(m+n.) 2 



(44) 

The exact energy splitting is obtained by numeri- 
cal diagonalization, in quadruple precision, of the 
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-log(AB) 



e , = 0.085: Exact 
Landau 



0.4 0.6 
E 



Figure 1: Negative logarithm of AE vs mean en- 
ergy E of pairs of almost degenerate consecutive 
energy levels. Inverted harmonic oscillator with 
h eff = 0.085 and 6= 1.5. 




Figure 2: Tunneling amplitude A vs mean energy 
E of pairs of almost degenerate consecutive energy 
levels. Inverted harmonic oscillator with h e ff = 
0.085 and b = 1.5. 



quantum Hamiltonian. We took a 4800 x 4800 ma- 
trix, thereby achieving 31 valid digits for the lower 
levels that we consider. In figure 1 we plot the neg- 
ative logarithm of AE function of the mean 
energy E of pairs of almost degenerate consecu- 
tive energy levels. We note a very good agreement 
between the exact and the Landau splittings. To 
resolve the differences, in figure 2 we show the tun- 
neling amplitude A (= the expression (29) without 
the exponential tunneling factor) as a function of 
the mean energy E of pairs of almost degenerate 
consecutive energy levels. 

6 Quartic potential 

In this section we consider the well known quartic 
potential, given by 

-2Bx 2 + Ax A , 



V(x) 



(45) 



where the parameters A and B are related to the 
potential barrier Vo and to the position of the min- 
imum xq by 



B 2 



x o = \ -r 



By using the following reduced variables 

x - E h 
— ' E = W 



x 



h, 



eff 



Xq 







(46) 



(47) 



the quantum Hamiltonian operator of the system 
can be written 



H = - 



2 Ox 2 



2x 2 + x A . 



(48) 



Let be a state of the quartic potential. It can 
be written in terms of the eigenstates \n) of the 
harmonic oscillator 



IV') = c "l n ) 

n=0 



(49) 



The Schrodinger equation for the harmonic oscil- 
lator is 

- f l^F + T* 2 I"> = EM ■ <50) 

By introducing the creation and annihilation op- 
erators 



1 



1 



(ujx — ip) 



(ujx + ip), 

which have the following properties 



(51) 



a + \n) = \Jn + l|n + 1) , 
a\n) = \/n\n — 1) , 



(52) 
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we get the matrix elements of the quantum Hamil- 
tonian of the quartic potential 



Hn 



= (n\H\m) 

+Km (^(2n 2 + 2n+l) 

■ in 



(53) 



x -h, 



ef 



OJ 1 

T + - 

4 uj 



+ ^(2n + 3) 



+0n+4,m I ^ 



:v /(n + 4)(n + 3)(n + 2)(n + 1) 



We calculate numerically, in quadruple precision 
(32 decimal digits), the energy levels of the sys- 
tem in the basis of the harmonic oscillator. For 
numerical purposes we took oj = 2 and the dimen- 
sionality 1000 x 1000, thereby again achieving 31 
valid digits for the lower levels that we consider. 

In this way we can compare the exact energy 
splittings with the semiclassical ones, which are 
obtained by using our splitting formula (29). For 
the quartic potential we have 



*o(0) 



V2 / ^/v{x)-Edx 



2^2 



er 2 (0) = - 



b 2 + E 







F(k) - bE(k) 



V" dx 



24V2dEJ ^V{x) - E 



24V2a 2 b 



(2 + 3E)b 2 + E 

VTTEb 4 



F(k) 



(54) 



^4(0) = - 



1152>/2 



d 2 



f 

Jo 



V"" dx 



dE 2 J y/v(x) - E 
7 d 3 r a V" 2 dx 



48\/2 



5 dE 3 Jo ^/V{x) - E 
[ZiF(k) + Z 2 E(k)] , 



where 



a = ^1-Vl + E, b = y/l + Vl + ti, k=( ^I 



Moreover we get 



rb I 
J a P 

f b °3 
Ja P 



dx 



dx = 



1 F(m) 
1 d 2 f b V" 



48V2dE 2 J a y/E-V 
1 



dx 



12^2 



f 

J a 



P 



dx = 



2304^ 



[^iF(m) + K 2 E(m)] , 
(55) 

7 d 4 r b V" 2 dx 

bdE^Ja y/E~^V 



d 3 



f 

J a 



V"" dx 



1440^ 



dE 3 J a VE~^V 
[QiF(m) + Q 2 E(m)] , 



where m 



b' 2 -a' 2 



and the functions F(m) and 



E(m) are the complete elliptic integrals of 1st and 
2nd kind, respectively, defined by (Abramowitz 
and Stegun 1972) 



dx 



^(l-x 2 )(l-mx 2 ) 



E(m) 



rra' 



dx . 



(56) 



(57) 



'0 \/l — x 

The functions Zi, Z 2 , Ki, K 2 , Q\ and Q 2 are de- 
fined as follows: 
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z 2 = 



K, = 



Qi 



Q2 



480(1 + E) 3 a% 5 
x (56b 2 + 135Eb 2 + 8WE 2 b 2 - lQ5E 3 b 2 

+112 + 333£ - 1050£ 2 - 375£ 3 ) , 

224 + 603E + 570£ 2 - 705£ 3 
480(1 + E) 3 a% 5 ' 

(1 - 7E) 



16(1 + E) 2 a 2 b 3 ' 

4 + 5E + 9E 2 
16(1 + E) 2 a A b 3 ' 



(58) 



128(1 + E) 9 / 2 a% 9 
x (b 2 (672 + 2263£ + 2730£ 2 + 9595£ 3 

-2008£ 4 - 1136S 5 - 528£ 6 ) 

+ (728-E + 2454.E 2 + 3755£ 3 

+ 11008S 4 - 349.B 5 + 528£ 6 )) , 

3(56 + 209£ + 289.E 2 + 195.E 3 - 165£ 4 ) 
8(1 +EYa% 7 



In most of the above manipulations we have used 
the Mathematica software. Please note that the 
quantities (13) needed in (29) are strongly diver- 
gent at the turning points, but all the expressions 
in (29) can be made finite by taking partial deriva- 
tives with respect to E of certain finite expres- 
sions as in (54) and (55), implying that in our for- 
mula (25) then all quantities are convergent and 
finite (c.f. Bender et al 1977, Robnik and Salas- 
nich 1997a,b). 

In figure 3 we show the tunneling amplitude A 
as a function of the mean energy E. We compare 



0.02s 



0.026 



0.024 



A 0.022 



h e{ = 0.03; Exact + 
Landau A 
Semiclassical up to fr^f * 



A * 
A 



* 



0.014 

-1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 
E 

Figure 3: Tunneling amplitude A vs mean energy 
E of Quartic potential with h e ff = 0.03. 



the exact results with the semiclassical ones at 1st 
(Landau), 3rd and 5th order in H. We observe 
that, as expected, there is a better agreement by 
increasing the energy because at high energy the 
classical momentum of the particle is large and 
thus the de Broglie wavelength sufficiently small 
for the semiclassical methods to be applicable. In 
figure 4 (5) we plot the tunneling amplitude A as 
a function H e f f for the first (fourth) pair of almost 
degenerate consecutive energy levels. As shown 
also in table 1, the semiclassical results approach 
the exact ones by increasing the perturbative order 
in h. Note that at the 5th order in h the agreement 
with exact result is up to the 8th digit. 
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A 0.045 



0.035 
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3 



Scmiclassical up to h e j * 
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* 
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0.025^ 

0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 



Figure 4: Tunneling amplitude A vs h e ff for the 
first pair of almost degenerate consecutive energy 
levels. Quartic potential. 



Figure 5: Tunneling amplitude A vs h e ff for the 
forth pair of almost degenerate consecutive energy 
levels. Quartic potential. 
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Table 1: Energy splitting AE for almost degenerate pairs at fi e ff = 0.03. E s is the exact mean energy, 
AE Landau, AE h 3 an d AE^s are the semiclassical results at 1st (Landau), 3rd and 5th order, respectively. 
In the last column of the upper table: exact energy splitting in units of mean level spacing. 



Pair number 


E s 




AE n /(El l+1 - E 1 }) 


1 


-0.9578013510838623 


9.232(2)E-28 


1.291(6)E-26 


2 


-0.8743363182686136 


6.533698E-25 


9.140541E-24 


3 


-0.7923227236203907 


2.197647599E-22 


3.074474586E-21 


4 


-0.7118393444104010 


4.66709587543E-20 


6.52919406381E-19 


5 


-0.6329771419061873 


7.0097983270055E-18 


9.8065981172742E-17 


6 


-0.5558426330893027 


7.901944841994295E-16 


1.10546971247045E-14 


7 


-0.4805627052945287 


6.920784581323967E-14 


9.68206928062435E-13 


8 


-0.4072917698084478 


4.806630946488390E-12 


6.72440144371474E-11 


9 


-0.3362229365299848 


2.675674688376366E-10 


3.74322699989559E-09 


10 


-0.2676066496781732 


1.196767019706533E-08 


1.67425832453064E-07 


11 


-0.2017846855780184 


4.273586943459416E-07 


5.97868122857108E-06 


12 


-0.1392610433553887 


1.195177462312803E-05 


1.67203455862220E-04 


13 


-8.088844930248588E-02 


2.500127371158503E-04 


3.49763905139656E-03 


14 


-2.855578436131249E-02 


3.352190557610379E-03 


4.68966210972998E-02 



^E Landau (E s ) 


AE h s(E s ) 


AE n5 (E s ) 


8 


.56386072023299E- 


28 


9 


33424576483334E- 


28 


9 


16991849898935E-28 


6 


34511775769322E- 


25 


6 


53796204696428E- 


25 


6 


53332080224002E-25 


2 


15758453279967E- 


22 


2 


19798578660872E- 


22 


2 


19763766003663E-22 


4 


60423712282166E- 


20 


4 


66737049838592E- 


20 


4 


66709238092202E-20 


6 


93436277656352E- 


18 


7 


01000176474235E- 


18 


7 


00979708264343E-18 


7 


83077639810467E- 


16 


7 


90207840692513E- 


16 


7 


90194445833250E-16 


6 


86703907509728E- 


14 


6 


92086233123417E- 


14 


6 


92078451430764E-14 


4 


77381431395796E- 


12 


4 


80667193943203E- 


12 


4 


80663098878842E-12 


2 


65943820818162E- 


10 


2 


67569528644699E- 


10 


2 


67567477543535E-10 


1 


19031027734624E- 


08 


1 


19677788146677E- 


08 


1 


19676714995564E-08 


4 


25353592752249E- 


07 


4 


27365612082292E- 


07 


4 


27358920983707E-07 


1 


19079278916304E- 


05 


1 


19524083756918E- 


05 


1 


19518371126262E-05 


2 


49835207740361E- 


04 


2 


50142635772383E- 


04 


2 


50058564705645E-04 


3 


46270869995442E- 


03 


3 


39489990444765E- 


03 


3 


33278741702504E-03 
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7 Conclusions 

In this work we have taken the first step towards 
a systematic improvement of the Landau formula 
(Landau and Lifshitz 1997), which is the semiclas- 
sical leading order energy level splitting formula 
for pairs of almost degenerate levels in double well 
potentials. We have developed the algorithm for 
the semiclassical fi expansion series to all orders 
for the tunneling amplitude A (of equation (1)), 
and thus calculated explicitly the quantum correc- 
tions up to the 5th order. We have compared the 
semiclassical predictions with the exact results ob- 
tained numerically, in quadruple precision in case 
of the quartic double well potential. Our approach 
is based on the usual WKB expansion in one- 
dimensional potentials. Thus the calculation of 
higher corrections can in principle be continued 
by the same method, although the structure of 
higher terms increases in complexity very quickly. 
We have also shown what happens in cases where 
the assumption implicit in the Landau formula 
(namely the linearity of the potential around the 
turning points) is not satisfied: We get a different 
result even in the leading semiclassical order, and 
this has been shown for the double square well po- 
tential. We should stress that the Landau formula 
(1997) is indeed quite good approximation since it 
always yields the correct order of magnitude (the 
exponential tunneling factor is always exact) and 
even the tunneling amplitude is correct within the 
5-50 % . 

It is our goal to work out a more direct WKB ap- 
proach to the solution of the multi-minima prob- 
lem, by the contour integration technique, based 
on requiring the single valuedness of the eigenfunc- 
tion, as has been done in the case of a single well 
potential in (Robnik and Salasnich 1997a, b). This 
is our future project. 

Finally we should mention important applica- 
tions e.g. in the domain of molecular physics 
(Herzberg 1991, Landau and Lifshitz 1997, Cohen- 
Tannoudji et al 1993). For example, the NH3 
molecule can be described by a quasi one- 
dimensional potential V(x) as a function of the 
perpendicular distance x of the nitrogen N atom 



from the H% plane, and as such it has the dou- 
ble well form, with the top of the potential bar- 
rier at x = 0. In order to calculate the energy 
level splitting of the doublets of vibrational modes 
one needs exactly our theory. Another example 
is the torsional motion of C2H4 molecule, where 
again we encounter an effectively one-dimensional 
double well potential. In case of C 2 Hq molecule, 
we find three potential wells where the tunneling 
effects again determine the splittings of energy 
triplet levels and a generalization of our theory 
would give an improved estimate of the splittings. 

Acknowledgments 

This work was supported by the Ministry of Sci- 
ence and Technology of the Republic of Slovenia 
and by the Rector's Fund of the University of 
Maribor. L.S. was partially supported by the MI41 
project of the Italian INFN. 

Appendix 

The semiclassical splitting formula is obtained by 
inserting the semiclassical wavefunction ipo(x) for 
the classically forbidden region \x\ < a 



i>o{x) 



C 



exp 



dx -\- Oeven + O'odd 



O'even = " a 2fc , 
k=l 

°2k(x)= I" a' 2k (\p(0\)d^ 



(59) 



&odd = 11 a 2k+l , 
k=l 

d 2k+1 (x)= f a' k (\m\)dt, 

into the general expression (8). The potential 
V{x) can always be written in the way where the 
maximum of the barrier is at the point x = 0. 



Considering that, and inserting (59) into (8) we 
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get 
AE 



2hC 2 
m\p\ 

x exp 



(|P| + Kven + °odd) 

J \p\ dp + a even + a odd 



(60) 



The exponent of all terms except the leading one 
can be expanded in the fi power series. The ex- 
pression valid up to the 5th order in % expansion 
reads as 



AE = 



2hC 2 



m 



l + ft2<r 2 (0) 



4(0) 



+ 2af(0) + 2a 3 (0) 



+h 3 ( 2d A (0) + |a|(0) + 4a 2 (0)a 3 (0) 



4(0) 



'4(0) 



The same combination multiplied by &2 and a 2 
can be also found in the brackets after T? and fi 4 . 
In the bracket after fi 4 one more combination like 
that can be found 



TA^ {x)+2al{x)+ m +2az{x) m 



2a' 5 (x)a' (x) + 2a' 4 (xy i (x) 
a' (x) 

2a' 3 (x)a' 2 (x) + a'l(x) 



+ 

+2a 3 (x) 
= 0. 



a' (x) 

2a' 3 {x)a' (x) + 2a' 2 (x)a[(x) + ~a' 2 \x) 
a' (x) 



(62) 

In this way the simplified splitting formula valid 
up to the 5th order expansion is found 



+K 



4 (°i(o) , O 4(o) / s2/ 



(2a 2 (0)|f|| + 2a5(0) + ^ 2 4 (0) 
+ h 4 f 2a|(0) + 4<7 2 (0)a 4 (0) + 4a|(0)a 3 (0) 



x exp 



2 



dx 



Using the recursive relation for a' k we observe that 
certain combinations in the tunneling amplitude in 
the upper equation are the integrals of functions 
identical to 0, evaluated in x = 0. In the term 
that stands with ft 2 such a combination is 



d fa' 2 (x) 



+ 2a 3 (x)^j = 



dx \a' (x) 

= 2a' 3 (xy (x) + 2a 2 {x)a[{x) + a' 2 '{x) 
a' (x) 



--0. 
(61) 



AE 



l + h2a 2 (0) +h 2 2dl(0) 

+S 3( 2i4(0) + l 5 3 (0) + |0|) 



(^^(0) + 4a 2 (0)a 4 (0) 



x exp 



2 /" a 



(63) 



and by taking into account the even symmetry of 
the potential V{x) due to which the first deriva- 
tives of all odd <Tfc evaluated at x = are zero, we 
obtain equation (24). 
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